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A statistical analysis of the elastic unpolarized electron proton scattering data shows that, at 
large momentum transfer, the size and the e dependence of the radiative corrections, as traditionally 
calculated and applied, may induce large correlations of the parameters of the Rosenbluth fit, which 
prevent a correct extraction of the electric proton form factor. Using the electron QED structure 
(radiation) function approach the cross section of elastic electron-proton scattering in leading and 
next-to leading approximations is calculated and expressed as a correction to the Born cross section, 
which is different for the electric and the magnetic contribution. When properly applied to the data, 
it may give the solution to the problem of the discrepancy of the polarized and unpolarized results 
on electron proton scattering. 

PACS numbers: 



I. INTRODUCTION 



The experimental determination of the elastic proton electromagnetic form factors (FFs) at large momentum transfer 
is presently of large interest, due to the availability of electron beams in the GeV range with high intensity and high 
polarization, large acceptance spectrometers, hadron polarized targets, and hadron polarimeters. The possibility of 
extending the measurements of such fundamental quantities, which contain dynamical information on the nucleon 
structure, has inspired experimental programs at JLab, Frascati and at future machines, such as GSI, both in the 
space-like and in the time-like regions. 

The traditional way to measure proton electromagnetic FFs consists in the determination of the e dependence of 
the reduced elastic differential cross section, which may be written, assuming that the interaction occurs through the 
exchange of one-photon, as [l|: 

af e r(<?,Q 2 ) = e(l + T) 



e= [1 + 2(1 + T)tan 2 (6»/2)]- 1 , 

where a — 1/137, r — Q 2 /(4to 2 ), Q 2 is the momentum transfer squared, to is the proton mass, E and 9 are the 
incident electron energy and the scattering angle of the outgoing electron, respectively, and Gm(Q 2 ) and Ge(Q 2 ) are 
the magnetic and the electric proton FFs and are functions of Q 2 , only. Measurements of the elastic differential cross 
section at different angles for a fixed value of Q 2 allow Ge(Q 2 ) and Gm{Q 2 ) to be determined as the slope and the 
intercept, respectively, from the linear e dependence IjTJl. 

High precision data on the ratio of the electric to magnetic proton FFs at large Q 2 have been recently obtained Q 
through the polarization transfer method Q. Such data revealed a surprising trend, which deviates from the expected 
scaling behavior previously obtained through the measurement of the clastic ep cross section according the Rosenbluth 
separation method 4]. New precise measurements of the unpolarized elastic ep cross section |5j and re-analysis of 
the old data 0,0 confirm that the behaviour of the measured ratio R(Q 2 ) — [iGe(Q 2 )/Gm{Q 2 ) = 2.79 is the 
magnetic moment of the proton) is different depending on the method used: 

• Scaling behavior for unpolarized cross section measurements: R(Q 2 ) — 1; Gm(Q 2 ) has been extracted up to 
Q 2 ~ 31 GeV 2 [8j and is often approximated, for practical purposes, according to a dipole form: Gr>{Q 2 ) = 
(1 + Q 2 /0.71 GeV 2 )- 2 ; 

• a strong monotonical decrease from polarization transfer measurements. 

R(Q 2 ) = 1 - (0.130 ± 0.005){Q 2 [GeV 2 ] - (0.04 ± 0.09)}. (2) 

The ratio deviates from unity as Q 2 increases, reaching a value of ~ 0.35 at Q 2 = 5 GeV 2 0. 

This puzzle has given rise to many speculations and different interpretations jjj llOl lll| , suggesting further experi- 
ments. In particular, it has been suggested that the presence of 2"/ exchan ge c ould solve this discrepancy through its 
interference with the main mechanism (I7 exchange). In a previous paper [12J it was shown that the present data do 
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not give any evidence of the presence of the 27 mechanism, in the limit of the experimental errors. The main reason 
is that, if one takes into account C-invariance and crossing symmetry, the 27 mechanism introduces a very specific 
non linear e dependence of the reduced cross section |l3l ITU. ll "J| , whereas the data do not show any deviation from 
linearity. 

Before analyzing the data in a different perspective, we stress the following points: 

• No experimental bias has been found in both types of measurements, the experimental observables being the 
differential cross section on the one hand, and the polarization of the outgoing proton in the scattering plane 
(more precisely the ratio between the longitudinal and the transverse polarization), on the other hand. 

• The discrepancy is not at the level of these observables: it has been shown that constraining the ratio R 
from polarization measurements and extracting Gm(Q 2 ) from the measured cross section "the magnetic FF is 
systematically 1.5-3% larger than had been extracted in previous analysis" , inside the error bars |l6(. 

• The inconsistency arises at the level of the slope of the e dependence of the reduced cross section, which is 
directly related to Ge(Q 2 ), i.e. the derivative of the differential cross section, with respect to e. The difference 
of such slope, derived from the two methods above, appears particularly in the latest precise data One 
should note that the discrepancy appears in the ratio R, whereas Gm{Q 2 ) decreases more than one order of 
magnitude from Q 2 =l to 5 GeV 2 . 

Radiative corrections to the unpolarized cross section can reach 30-40% at large Q 2 . RC are calculated as a global 
factor which is applied to the number of detected elastic events, a meas . As a rule, they depend on the kinematical 
variables, as e and Q 2 . They are traditionally applied to the unpolarized differential cross section, following a 
prescription which includes only leading order contributions [TtI |2|| : 



W 4) ; r.u 



4 V m 



where f(6) is function only of the scatte ring angle 8. 

As noted in the original papers [H Hf, when AE 0, <7 meas becomes negatively infinite, whereas physical 
arguments require that it should vanish. The authors stated already that this problem would be overcome taking into 
account higher order radiative corrections. 

In recent experiments E is large and the experimental resolution is very good (allowing to reduce AE). Moreover, 
multiple photon emission from the initial electron, shifts the momentum transferred to the proton to lower values, 
increasing the cross section. Therefore S becomes sizeable and one can not safely neglect higher order corrections. 

A complete calculation of radiative corrections should take into account consistently all different terms which 
contribute at all orders (including the two photon exchange contribution) and their interference. 

However, several approximations are made, which may not be safely extrapolated to the conditions of the present 
esperiments. In particular in the calculation of Ref. |17| . the consideration of hard collinear photon emission (where 
the radiative photon is emitted along the direction of the incident or outgoing electron) is not complete. Moreover 
higher order RC, pair production as well as vacuum polarization are not included. 

We calculate here the cross section of elastic electron-proton scattering in leading and next-to leading approximation 
using the electron QED structure (radiation) function approach, which takes into account any number of real and 
virtual photons, emitted in collinear kinematics, at all orders in QED. This approach was previously applied to 
unpolarized e + e~ scattering ^8|, to deep inelastic scattering 0] and, more recently, to polarization observables in ep 
elastic scattering [53, Hi). It was found that the correction is lower than 1% to polarization observables, as expected. 
However, the effect on the polarized cross section was not investigated, and, in particular the effect on the slope of 
the reduced cross section as a function of e, which is the relevant quantity here. 

The purpose of this paper is to re-analyze the unpolarized data, with particular attention to the applied RC. 

The paper is organized as follows. In Section II we show that a large correlation exists between the two parameters 
extracted from the Rosenbluth fit at large Q 2 and analyze the existing data in this respect. A probable source of 
these correlations being found in the standard procedure taken for RC, we calculate RC for the cross section of elastic 
electron-proton scattering in frame of the structure functions (SF) approach (Section III). In Section IV numerical 
results are presented, which show that the correction to the measured cross section is different for the electric and the 
magnetic contribution, and therefore affects the slope of the reduced cross section and, in particular, the extraction 
of G E (Q 2 )- 



3 



2.5 



1 .5 
1 

0.5 



a 





1.1 0.8 

. 2 / 2a 2 



1 1.1 1.2 



2^ 2 



FIG. 1: Dependence of G 2 E /G 2 D versus G 2 M / fj, 2 G 2 D : (a) for Q 2 > 2 GeV 2 from Refs. p (triangles), @ (stars) and O (squares); 
(b) for Q 2 < 2 GeV 2 from Refs. |H (circles), and (squares). 



II. STATISTICAL ANALYSIS OF THE PRESENT DATA 



The starting point of this work is the observation of a correlation, which appears in the published FFs data extracted 
with the Rosenbluth method: the larger is G E , the smaller GL The dependence of G 2 E /G 2 D versus G 2 M /ii 2 G 2 D is 
shown in Fig. ^ for three recent data sets, at Q 2 > 2 GeV 2 PIlI ■ I n Fig- ffl 3 tw0 data sets at low Q 2 (Q 2 < 2 
GeV 2 ) are shown Whereas at low Q 2 , G 2 E /G 2 D seems constant and quite independent from G 2 M / ii 2 G 2 D , at 

large Q 2 an evident correlation appears. This is especially visible in the most recent and precise experiment, at large 
Q 2 where a linear fit of the ratio R{Q 2 ) as a function of Q 2 gives: 

R{Q 2 ) = (0.13 ± 0.11){Q 2 [GeV 2 ] + (0.57 ± 0.32)}. (4) 

where Q 2 is expressed in GeV 2 . 

Polarization data also showed a linearity of the ratio R, with the same slope (in absolute value), Eq. but with 
opposite sign. In this case, the ratio is measured directly, whereas according to the Rosenbluth method one extracts 
two (independent) parameters from a linear fit. A correlation between the two parameters could be induced by the 
procedure itself or could be a physical effect and have a dynamical origin. In the latter case, it should not depend on 
the experiment. 

It is known that at large Q 2 the contribution of the electric term to the cross section becomes very small, as the 
magnetic part is amplified by the kinematical factor r. This is illustrated in Fig. |2 where the ratio of the electric part, 
Fe = eG|;((5 2 ), to the reduced cross section is shown as a function of Q 2 . The different curves correspond to different 
values of e, assuming FFs scaling (thin lines) or in the hypothesis of the linear dependence of Eq. (thick lines). In 
the second case, one can see that, for example, for e = 0.2 the electric contribution becomes lower than 3% starting 
from 2 GeV 2 . This number should be compared with the absolute uncertainty of the cross section measurement. When 
this contribution is larger or is of the same order, the sensitivity of the measurement to the electric term is lost and 
the extraction of Ge{Q 2 ) becomes meaningless. Since the first measurements [2^, electromagnetic probes have been 
traditionally preferred to hadronic beams, as the electromagnetic interaction is exactly calculable in QED, and one 
can safely extract the information from the hadronic vertex. However, one has to introduce the radiative corrections, 
which become very large as the momentum transfer squared, Q 2 , increases. Radiative corrections were first calculated 
by Schwinger 26] and are important for the discussion of the experimental determination of the differential cross 
section. 
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FIG. 2: Contribution of the Ge(Q 2 ) dependent term to the reduced cross section (in percent) for e = 0.2 (solid line), e = 0.5 
(dashed line), e = 0.8 (dash-dotted line), in the hypothesis of FF scaling (thin lines) or following Eq. @ (thick lines). 



The measured elastic cross section is corrected by a global factor Cr, according to the prescription |17| : 

^Born r i ^meas /r\ 

a red - ^R^red • l^J 

The factor Cr contains a large e dependence and a smooth Q 2 dependence, and it is common to the electric and 
magnetic parts. At the largest Q 2 considered here this factor can reach 30-40%, getting larger when the resolution is 
higher. If one made a linear approximation for the uncorrected data, one might even find a negative slope starting 
from Q 2 > 3 GeV 2 0. 

In Fig. |3| we show the Cr dependence on e for different Q 2 and from different sets of data. One can see that Cr 
increases with e, rising very fast as e — > 1. It may be different in different experiments because its calculation requires 
an integration over the experimental acceptance. 

The Rosenbluth separation consists of a linear fit to the reduced cross section at fixed Q 2 , where the two parameters 
are G 2 E and G 2 M . The multiplication by a factor which is common to the electric and magnetic terms, see Eqs. (11151 . 
and depends strongly on e, induces a correlation between these two parameters. In order to determine quantitatively 
how large this correlation is, we have built the error matrix for the Rosenbluth fits to the different sets of data available 
in the literature. 

At fixed Q 2 the reduced cross section, normalized to G 2 D , has been parametrized by a linear e dependence: 
c^"J n /G 2 D = ae + b. The two parameters, a and b, have been determined for each set of data as well as their 
errors a a , and the covariance, cov(a, b). The correlation coefficient £ is defined as £ = cov{a, b)/a a cTb and is shown 
in Fig. Has a function of the average of the radiative correction factor < Cr >, weighted over e. 

As the radiative corrections become larger, the correlation between the two parameters also becomes larger, reaching 
values near its maximum (in absolute value). Full correlation means that the two parameters are related through a 
constraint, i.e. it is possible to find a one-parameter description of the data. The data shown in Fig. 01 correspond to 



5 



2£ 



1 .6 



.5 



4 



1 .5 - 



1 .2 




1.1 







0.2 



0.4 



0.6 



FIG. 3: Radiative correction factor applied to the data at Q 2 =3 GeV 2 (squares) from Ref. [Zj, at Q 2 —4 GeV 2 (triangles) and 
5 GeV 2 (reversed triangles) from Ref. Q, and at Q 2 =0.32 GeV 2 from Ref. [24| (circles). The lines are drawn to guide the eye. 



those sets of experiments where the necessary information on the radiative corrections is available. The correlation 
coefficient itself can be calculated for a larger number of data and it is reported in Table [I] 
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Q 2 [GeV] 




Ref. 


Q 2 [GeV] 




Ref. 


2.6400 


-0.8823 


M 


0.2717 


-0.7258 


[24] 


3.2000 


-0.8973 




0.2911 


-0.7818 




4.1000 


-0.9060 




0.3105 


-0.7085 




1.7500 


-0.8693 


\£ 


0.3493 


-0.7683 




2.5000 


-0.9141 




0.3881 


-0.7417 




3.2500 


-0.9242 




0.4269 


-0.7093 




4.0000 


-0.9178 




0.4657 


-0.7381 




5.0000 


-0.8940 




0.5045 


-0.8126 




1.0000 


-0.9918 


[22] 


0.5433 


-0.7646 




2.0030 


-0.9915 




0.5821 


-0.8076 




2.4970 


-0.9910 




0.6209 


-0.8061 




3.0070 


-0.9878 




0.6598 


-0.8137 




0.1552 


-0.6761 


[24J 


0.6986 


-0.8713 




0.1785 


-0.6788 




0.7374 


-0.8145 




0.1940 


-0.6915 




0.7762 


-0.8512 




0.2329 


-0.7177 




0.8538 


-0.7612 





TABLE I: Correlation coefficient £ = cov(a,b)/a a ai 1 for different sets of data. 



At low Q 2 a correlation still exists, but it is smaller. For the data from Ref. [2^ the radiative corrections are of 
the order of 15%, seldom exceed 25% and correspond to e < 0.8. This allows a safer extraction of the FFs. 

Fig. 01 shows that, for each Q 2 , the extraction of FFs by a two parameter fit may be biased by the e dependence 
induced by the radiative corrections. Whatever the precision of the individual measurements is, the slope of the 
reduced cross section is not sensitive to Ge(Q 2 ) at large Q 2 . The Q 2 dependence is therefore driven by Gm(Q 2 ), 
which follows a dipole form. For each Q 2 a nonzero value of the ratio R will lead to an apparent dipole dependence of 
Ge(Q 2 )- Therefore experiments based on this method, will always give a Q 2 dependence of Ge(Q 2 ) which is driven 
by Gm(Q 2 ), i-e. follow approximately a dipole behavior. 

III. CALCULATION OF RADIATIVE CORRECTIONS 

It is known that the process of emission of hard photons by initial and scattered electrons plays a crucial role, 
which results in the presence of the radiative tail in the distribution on the scattered electron energy. We give here 
a simplified example of how a different calculation of radiative corrections can affect the electric and magnetic part 
of the unpolarized cross section, and change, in particular, its e dependence. The main point of interest here is to 
show the very sharp dependence of the initial state emission on the inelastic tail of the scattered electron energy 
spectrum. A more extended version of this calculation and its application to polarization observables, including two 
photon exchange, is given elsewhere |28j . 

The aim of this paper is to drive the attention to the sensitivity of the electric FF to the procedure used for its 
extraction from the data, and to focus the attention on how radiative corrections are applied on the unpolarized cross 
section. 

The structure functions (SF) approach extends the traditional one |17||. taking precisely into account the contri- 
butions of higher orders of perturbation theory and the role of initial state photon emission. The cross section is 
expressed in terms of SF of the initial electron and of the fragmentation function of the scattered electron energy 
fraction. Experimentally the detection of the scattered electron does not allow to separate the collinear photon emis- 
sion. Therefore, one integrates in a range of the scattered electron energy. This is equivalent to set the fragmentation 
function to unity, due the well known properties of this formalism 19] . 

It is known that initial state photon emission is more important than final state photon emission, due to the effect 
of decreasing Q 2 . Proton emission (real and virtual) is essentially smaller than the electron ones, and can be included 
as a general normalization. Vacuum polarization, which has been often neglected in previous analysis, here is taken 
into account. 
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FIG. 4: Correlation coefficient, f, as a function of the radiative correction factor < Cr >, averaged over e, for different sets of 
data: from Ref. [2^1 (circles), from Ref. 3 (triangles) and from Ref. (squares). 



The four momentum transfer squared can be written as: 

2 2£ ,2 (l-cos( 



where p is the recoil factor: 
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1 



In an experiment, the selection of elastic scattering requires the integration of the events in the elastic peak, and 
the rejection of inelastic events. We parametrize the cut on the energy of the final electron E' , selecting events with 
E' > c/p, where c is the 'inelasticity' cut, c < 1 (for the present numerical application we choose c = 0.97). 
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Due to the properties of SF method, radiative corrections can be written in form of initial and final state emission, 
although gauge invariance is conserved. This form obeys the Lee-Nauenberg-Kinoshita theorem, about the cancel- 
lation of mass singularities, when integrating on the the final energy fraction. This results in omitting the final 
(fragmentation) SF, i.e., in replacing the structure function associated with the final electron emission by unity. 

Therefore, the differential cross section, calculated in frame of the SF method, dc ^ n , can be written as |27j : 



da 



SF 



a 2 cos 2 (6>/2) 

~dVT ~ AE 2 sin 4 (6>/2) 



dzD(z) 



[i-n(Q= 



(i + H 

V 7T 



K 



(6) 



where K is an e-independent quantity of the order of unity, which includes all the non-leading terms, as two photon 
exchange and soft photon emission. More precisely the interference between the two virtual photon exchange amplitude 
and the Born amplitude and the relevant part of the soft photon emission i.e., the interference between the electron and 
proton soft photon emission, are included in the term K. This effect is not enhanced by large logarithm (charcteristic of 
SF) and can be included in non-leading contributions. The factor 1 + —K can be considered as a general normalization. 
It is calculated in detail in Ref. |28| . 

Here we focuss on the e-dependence of the differential cross section. The SF calculation, Eq. 0J, can be expressed 
in the form of a correction the Born reduced cross section (we omit RC of higher order) : 



r SF 



°Ved U + A ) 



(7) 



with 



A»-2{§('-!>-i(*-i> 



2 In 



l-z 



zo 



1 



p(l + r)(Z-l) 



1 {l + z 2 )dz 



1 - z 



(1) 



[i-n(z)]2 [i-n(i])2 



, L = In 



Q 2 



m e is the electron mass. The structure (radiation) function D(z) is 



D(z) = 



+ 0(f), (3 



2a 



Q 2 
ln^-1 

■mi 



The lower limit of integration, z , is related to the 'inelasticity' cut, c, necessary to select the elastic data: 

c 



p-c(p- 1) 



(8) 



(9) 



(10) 



The transfer momentum and recoil factor of the scattered electron after the collinear photon emission are, respectively, 
Q z and p z : 



Ql = 2E 2 z 2 {\ ~ cos9)/ Pz ; p z = 1 + z— (1 - cosi 

The kinematically corrected Born cross section for the scattered electron, <fi(z), is: 

1 



with 



e z z 2 p z (l + t z ) 

_ Ql i 



Ured{z), Vred(z) = T Z G M (Q Z ) + C Z G E {Q Z ) 



l + 2(l + r z )tan 2 (6>/2). 



(11) 



(12) 



(13) 



The vacuum polarization for a virtual photon with momentum q, q 2 = —Q 2 < 0, is included as a factor 1/[1 — n(Q )]. 
The main contribution to this term arises from the polarization of electron-positron vacuum: 



n(Q 2 ) 



3tt 



L 



(14) 
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The calculation requires a specific procedure for the integration of the SF D(z), which has a singularity at the 
upper limit of integration, Eq. © . 

The dependence of SF reduced cross section, Eqs. (|T114|1 . on e is shown in Figs. EJ'jbjC, for different values of 
Q 2 =l, 3, 5 GeV 2 , (solid lines). For comparison, the corresponding Born reduced cross section assuming also FFs 
parametrized in dipole form is shown as a dashed line, and the Born cross section, with FFs parametrized according 
to polarization measurements, (Eq. J2J as a dash-dotted line. 

One can see that SF corrections affect the e dependence of the cross section. Such effect is more important as 
Q 2 increases and for large e values. The relative difference of the SF reduced cross section with respect to the Born 
reduced cross section (both assuming dipole FFs), |A SF |, is shown in Fig. [SJi. For large values of e, the calculated 
reduced cross section can differ from the Born one by more than 7%, for c = 0.97. As both calculations assume dipole 
FFs, the source of the difference has to be attributed to how radiative corrections are calculated and applied. 

Let us stress that the main effect of this correction is to modify and lower the slope of the reduced cross sec- 
tion. This effect brings into qualitative and quantitative agreement FFs data issued from polarized and unpolarizcd 
measurements, as one can see from the comparison of the solid and dash-dotted lines in Figs. Et,b,c. 

Of course, the concrete value of the slope depends on the inelasticity cut. Taking 0.95 < c < 0.97, the slope given 
by the SF calculation is in complete agreement with the slope suggested by the polarization measurements. 

IV. CONCLUSIONS 

We reanalyzed the Rosenbluth data with particular attention to the radiative corrections applied to the measured 
cross section, and we showed from the (published) data themselves that at large Q 2 statistical correlations between 
the parameters of the Rosenbluth plot become so large that Ge(Q 2 ) can not be safely extracted. The method itself 
is biased at large momentum transfer because RC are applied as a global factor, which is the same for the electric 
and the magnetic contribution. Such factor contains a large e-dependence, which induces a strong correlation in the 
parameters of the linear e fit. 

Calculations of RC in frame of the SF method, which takes into account higher order of perturbation theory, 
show that RC from collinear hard photon emission affect the elastic ep cross section, in particular its e dependence. 
Similarly to the standard RC, they depend on the electron scattering angle and on the kinematical selection for the 
elastic events. On the opposite, they act differently on the electric and magnetic term of the cross section, changing 
the slope of the reduced cross section which is related to the electric FF. When applied to the polarized cross section, 
it has been shown that their effect is small on the relevant observables |2jJ, El ■ Therefore it is suggested here that 
such corrections, when properly applied to the experimental data, can bring into agreement the results on the proton 
FFs issued from unpolarized and polarized measurements. Moreover these corrections affect very little the linearity 
of the Rosenbluth fit, contrary to what is expected from two photon exchange |l3[ . 

A complete calculation should take into account consistently all different terms which contribute at all orders 
(including the two photon exchange contribution) and their interference |2S|. 

We confirm the conclusion of a previous paper Q which first suggested the polarization method for the determination 
of Ge(Q 2 ), due to the increased sensitivity of the cross section to the magnetic term at large Q 2 : 'Thus, there exist 
a number of polarization experiments which are more effective for determining the proton charge FF than is the 
measurement of the differential cross section for unpolarized particles' . 
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